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We present a general formalism based on the variational principle for finding the time-optimal 
quantum evolution of mixed states governed by a master equation, when the Hamiltonian and 
the Lindblad operators are subject to certain constraints. The problem reduces to solving first a 
fundamental equation (the quantum brachistochrone) for the Hamiltonian, which can be written 
down once the constraints are specified, and then solving the constraints and the master equation 
for the Lindblad and the density operators. As an application of our formalism, we study a simple 
one-qubit model where the optimal Lindblad operators control decoherence and can be simulated 
by a tunable coupling with an ancillary qubit. It is found that the evolution through mixed states 
can be more efficient than the unitary evolution between given pure states. We also discuss the 
mixed state evolution as a finite time unitary evolution of the system plus an environment followed 
by a single measurement. For the simplest choice of the constraints, the optimal duration time for 
the evolution is an exponentially decreasing function of the environment's degrees of freedom. 
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I. INTRODUCTION 

Recently, many works related to time optimal quantum 
computation have appeared in the literature [l|, 0, El, 0, 
H, H 0, B H, OH 03 ■ The minimization of physical time 
to achieve a given unitary transformation should provide 
a more physical description of the complexity of quan- 
tum algorithms. In a scries of previous works [l3| 
we established a general theory based on the variational 
principle to find the time optimal evolution between given 
initial and final pure states [l2j (paper I), and to find the 
time optimal way of generating a target unitary opera- 
tion for arbitrary initial states (l3j (paper II). In paper 
I we studied closed pure quantum systems driven by the 
Schrodinger equation and where the Hamiltonian is con- 
trollable within a certain available set. Paper II is an 
extension of paper I and is more relevant to subroutines 
in quantum computation, where the input may be un- 
known. The main result of our works is that, once the 
constraints for the Hamiltonian are given, one can sys- 
tematically derive a fundamental equation, the quantum 
brachistochrone, which can be always solved, at least nu- 
merically, for the time optimal Hamiltonian (without any 
further restricting assumptions, e.g. the adiabaticity of 
the quantum evolution). Here we extend our previous 
works and formulate a variational principle for the time 
optimal quantum control of open systems where the dy- 
namics is driven by a master equation in the Lindblad 
[II, EH form. 

Historically, quantum control theory of pure states has 
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been studied b y m any people (for a review of the sub- 
ject, see, e.g., [lal). Around twenty years ago, Peirce, 
Dahleh and Rabitz [ItJ considered a variational method 
to manufacture a wave packet as close as possible to a 
given target. Other authors (see, e.g., [IH and other ref- 
erences in [HI]) further investigated the variational meth- 
ods by optimizing the fidelity between the final state of 
the steered system and a given target state. The applica- 
tion to the optimal realization of unitary gates in closed 
systems was also studied (see, e.g., 0). For the mixed 
state case, the master equation in the Lindblad form has 
been used in [2(| (for other recent references see, e.g., 
[2ll . I22I , l23jj. However, while the main concern of these 
papers was the optimization of the fidelity or the purity 
of the quantum operations, here we focus the attention 
on the time optimality. 

The paper is organized as follows. In Sec. [II] we in- 
troduce the Markovian approximation and the master 
equation for an open quantum system and we discuss 
the related gauge degrees of freedom. In Sec. IIHI we set 
up the general variational formalism for the time optimal 
evolution of such quantum systems. The action consists 
of a time cost function to be minimized and of Lagrange 
multiplier terms which ensure the evolution under the 
master equation and certain constraints for the available 
Hamiltonian and Lindblad operators. We then derive the 
fundamental equations of motion. In Sec. lIVI we explicitly 
demonstrate our methods via the example of a one-qubit 
system by deriving the time optimal Hamiltonian, Lind- 
blad and density operators, which may represent cither 
a measurement or a decoherence process. In Sec. |V] 
we simulate the optimal operations derived in Sec. IIVI 
by the partial trace of a two-qubit self-interacting system 
with a controllable Hamiltonian and ancillary qubit. This 
corresponds to a repetition of short-time measurements. 
Sec. I VII instead, describes the mixed state evolution 
as a finite time unitary evolution of the system plus an 
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environment followed by a single measurement. Finally, 
Sec. lVIIl is devoted to the summary and discussion of our 
results. 



II. MASTER EQUATION 

In this paper, we address the problem of time optimal 
quantum control of open systems where the dynamics 
is described by a master equation in the Lindblad form 

p := C(p) = -i[H, p] + J2 UaPLl - \{L\L a , p}) (1) 



for the density operator p(t), where H(i) is the Hamil- 
tonian, L a (t) (a = 1,...,N 2 — 1) are the Lindblad op- 
erators, N is the dimension of the Hilbcrt space of the 
system and we use the notation A := dA/dt for time 
derivatives, [A, B] := AB — BA for the commutator and 
{A, B} := AB + BA for the anticommutator. The Hamil- 
tonian represents the unitary evolution part while the 
Lindblad operators express generalized measurements or 
decoherence processes due to the coupling of the system 
with an environment. The master equation is Markovian, 
i.e., the environment has no memory of the main physical 
system. It can be physically realized if the interaction be- 
tween the main system and the environment is weak and 
its typical time scale is small compared with that of the 
physical system 0, [H[ . A simple example discussed in 
Sec. V illustrates how the repetition of a short time uni- 
tary evolution and the partial trace (e.g., measurement) 
over the environment states (|39|) reproduces the master 
equation JT]). 

The evolution of p(t) is invariant under the following 
gauge transformations 



H » H + a(t) + I - PaWt], 



L a ^ L a + P a (t), 



and 



(2) 



(3) 



where a(t) is a real number, [3 a (t) are complex numbers, 
and the U a (,(t) form a unitary matrix with respect to its 
indices (i.e., J2 c u acU bc = J2 c U* a U cb = S ab ). The pa- 
rameter a(t) in ^ corresponds to the £7(1) gauge degree 



of freedom which was discussed in our first papers 12|, 13[ 
for the case of pure states. The gauge degrees of freedom 
Pa(t), instead, correspond to the fact that the operator 
H(t) in ([T|) is not just the free Hamiltonian of the re- 
duced system, but may contain coupling terms with the 
environment (see, e.g., equation in Sec. V). Finally, 
the U a t(t) represent the freedom of the choice of the basis 
for the HUbert space of the environment. 



From the gauge freedom ([2]), we can always choose 
H = H, L a = L a , where a tilde denotes the traceless 
part of an operator, A := A — (1/N) Tr A. Furthermore, 
by using the degrees of freedom of ((3]), we can also choose 
the L a to be mutually orthogonal. This is the gauge 
choice which we make throughout this paper. 



III. GENERAL FORMALISM 

Let us consider the problem of controlling a certain 
physical system governed by the master equation |T]) and 
of steering its transition between given initial and final 
quantum states in the shortest time. Mathematically 
this is a time optimality problem for the evolution of 
the density operator p(t) according to |T]) by controlling 
the Hamiltonian H{t) and the Lindblad operators L a (t). 

We assume that at least the 'magnitudes' of the Hamil- 
tonian and of the Lindblad operators are bounded. Phys- 
ically this corresponds to the fact that one can afford only 
a finite energy in the experiment, and that a finite level 
of decoherence is tolerated. Besides this normalization 
constraint, the available operations may be subject also 
to other constraints, which can represent either experi- 
mental requirements (e.g., the specifications of the ap- 
paratus in use) or theoretical conditions (e.g., allowing 
no operations involving three or more qubits). We then 
consider the following action for the dynamical variables 
p(t), H(t) and L a (t) 



S(p, H, L a , a, X a ) = J dt [L T + L M + L c ] , 



with 



L 



9p{(>,p) 



Y g p {L{p),L{p)Y 
L M :=Ti[a(p -£(/>))], 
Lc ■= X afa(H) + ^/i afe [Tr(4L 6 



(4) 



(5) 

(G) 
(7) 



where the traceless Hcrmitian operator a, the real func- 
tions A a and the complex functions p a b = M& a are -"- ,a " 
grange multipliers, and Tt p(t) = 1. As mentioned in 
the previous section, we assume (even before taking vari- 
ations of the action) that H and L a are traceless [26| . 
while the second term in Lc is a constraint which en- 
sures the normalization and the mutual orthogonality of 
the L a operators. 

The Lt term in the action (|2|) gives the time dura- 
tion for the evolution of p(t). The Riemannian metric g p 
above is assumed to belong to the family of monotone 
metrics on the space of density operators given by 



g p (A,B) := Tr[Ac(L,R)(B)], 



(8) 



where L and R are multiplication of p from the left 
and right, respectively, c{x,y) := \ / (yrj(x / y)) , and r) is 
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an operator- monotone function satisfying trjit^ 1 ) = r](t) 
[27l . |28| . The equations of motion derived below actually 
do not depend on the particular choice of g p . In fact, 
they are the same for all Lt which are constant on shell, 
i.e., when the master equation holds (and one may even 
choose Lt = 1). 

The Lc term in (0| generates the constraints 



and 



fa(H) = 



Tr(tf a L b ) = NjXb 



(9) 



(10) 



for the Hamiltonian and the Lindblad operators. The 
Lm term in @ guarantees that p(t) satisfies the master 
equation {l}. 

The other equations of motion are derived in the fol- 
lowing way. From the variation of S by p, with the help 
of {!]), we get the adjoint master equation 



([16]) is the same universal equation as for pure states [12| 
and for unitary operations flU ] . It can also be obtained 
by observing that the first two terms Lt and Lm in the 
action Q are invariant under an arbitrary infinitesimal 
non-Abelian transformation 



p^ p-i[x(t),p], 
a 1 > o- + i[x{t),a], 

H»H + X (t)-i[x(t),H}, 
L a ^ La- i[x(t),L a ], 



(17) 



where x(t) G su(iV). This is because (fT7j) does not 
change the master equation |T]) and the adjoint master 
equation (fTTj) . Therefore the variation of the entire ac- 
tion reduces to the variation of the constraint term Lc, 
which is easily checked to produce the quantum brachis- 
tochrone equation (fT6|) . The same derivation also holds 
for the case of pure states and unitary operations treated 
in our previous works [HI, [l3[ ■ 



(ii) 



for a' := a + c(L, R)(£(p))/ g p (C(p), C{p)) , where the su- 
peroperator £t j s defined by Tr[A^ C{B)] = Tr[tf (A)B] 
and its explicit form is tf(A) = i[H, A] + J2 a ( L lAL a - 
\{L\L a , A}). The variation of S by H and JT]) imply 



F = -i[p,a'}, 
where we have defined the operator 

dL c 



F(H) 



dH ' 



(12) 



(13) 



which will be important in the sequel. From the variation 
of 5 by L\, together with ([1]), we obtain 



traceless part of 



°~'L a p- -^L a {p, cr'} 



^ P-abLb- 
b 

(14) 



One can show that p a i> is actually diagonal, i.e. that 
Pab = PaSab ■ and therefore get the algebraic formula 



traceless part of 



<j'L a p- ^L a {p 7 cr'} 



= p a L a , (15) 



which is an eigenvalue equation with eigenvalues p a for 
the eigenvectors L a . Finally, combining (TTJ), pT|) . (fT2"|) 
and (|15p . we can eliminate the Lagrange multipliers to 
obtain the fundamental quantum brachistochrone equa- 
tion 



IF = [H, F] 



(16) 



We can thus obtain the time optimal H{t) and L a {t) 
separately. One first solves (fir?)) to obtain the optimal 
H{t) and then solves (1]), (JTTJ) , (JT2J) and (TT5D to nnd tnc 
optimal L a (t). The quantum brachistochrone equation 



IV. A ONE-QUBIT EXAMPLE 

Let us now discuss as an explicit example a one-qubit 
model where the Hamiltonian is subject only to the nor- 
malization constraint 

f (H):=±(TrH 2 -Nu 2 ) = 0, (18) 

where w is a given constant. In this case F = XqH and 
the quantum brachistochrone equation (JTHJ) becomes triv- 
ial, giving H = const and A = const (see Sec. Ill in 
[13|). Our problem then reduces to that of solving the 
master equation |T]) and the adjoint master equation pip 
together with the algebaric equations (fT2")) and (|15|) for 
p(t), L a (t) and tx'(i). In the Pauli basis {a x , a y , a z } these 
can be rewritten as equations for three-dimensional vec- 
tors. Namely, if we parametrize the states as 



p(t) := _[l + r (i)-<7], (19) 
a'(t) := s(t) ■ cr, (20) 
and the Hamiltonian and the Lindblad operators as 

H := h er, (21) 



L a it) := l a (t) ■ cr, 



(22) 



where r, s and h £ M 3 and /„ e C 3 , the master equation 
([1]) and the adjoint master equation (fTTj) can be rewritten 
as 

r = 2[hxr+Y J (M(* ' r ) 1 * )~ \ l ? r + il x **)]> ( 23 ) 
ie{i a } 

s = 2[/ixs-^ (Re((i • s)l*) - \l\ 2 s)}. (24) 
ie{i a } 



Moreover, the algebraic equations (jl2|) and ([15)1 read 

r x s = Xoh = const (25) 
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and 

K(r,a)l a =v a l a , (26) 
where K(r, s) is the self-adjoint matrix 

Kj k (r, s) := rjs k + r fc Sj - 2i ^ e jU si (27) 

Z— x,y,z 

and z/ a := 2(r • s + fi a ) is a real number. 

Because of the constraints (fT8"|) and (f^5|) . the compo- 
nents of the Hamiltonian h are given by 



4-1 . rx» 

wn if r x s = 0, 



if rxs^O, 



(28) 



where n is an arbitrary unit vector. The components of 
the Lindblad operators l a are determined as eigenvectors 
of the eigenvalue equation (j2"6")l with the constraints (flT))) . 
i.e. |Z | = 7 a . 

At a given instant, we parametrize r(t) and s(i) as 



r(t) 
s(t) 



r [ cos -e 3 + sin -ei 



s ( cos -e 3 - sin -ei 



(29) 
(30) 



where {ei(i), e2(i), e 3 (i)} are orthonormal vectors, so 
that r ■ s = rs cos 0, with r £ [0, 1] and 9 £ [0, 7r] . We can 
then rewrite the matrix if in (l27l) as 



of the Lindblad operators arc different. Note that in this 
particular case the Lq operator, which corresponds to a 
projective measurement along e^, is not effective for the 
state evolution while the amplitude damping L± play a 
significant role in (|35|) . The case h = is depicted in 
Fig. [T] In the case when rxs^O the coupled equations 




FIG. 1: Analytical, time optimal evolution of p(t) (arrows) in 
the Bloch sphere for the case of r x s = 0, 7+7^0, 7- = 
and p(0) = |T)(T|- Also shown (thick solid meridian curve) 
the optimal pure state evolution between the north and south 
poles 0. 



(K jk ) = 2s 



l COS : 







i sin ■ 



([2"3")l . and (|2^|) can be solved numerically. We depict 
(31) a family of optimal trajectories from a mixed state to a 
pure state in the Bloch sphere in Fig. 2. 



If the conserved vector satisfies r x s = 0, we can see 
6 = (i.e. r = re^), and the components of the Lindblad 
operators are given by (|3ip as the following vectors 



7± 

l±(t) = - 7 =(e 1 ±ie 2 ), 
lo(t) = 7 e 3 . 



(32) 
(33) 



To simplify the analysis, we can move to the interaction 
picture by the transformation p' = UqpUq := |(l+r' • 
<t) with Uo(t) = Tcxp(— i J H(t)dt), so that the master 
equation |23|) becomes 



-2[( 7 i+7-y-(7l-7^)e 3 ], 



(34) 



which guarantees e 3 = 0. Therefore we obtain the fol- 
lowing solution for the Bloch vector in the interaction 
picture 



r'(t) 



l 2 - 



7+ +7- 



r'(0) 



-2(7+ +7_)« 



e 3 . 



(35) 

If the magnitudes of the Lindblad operators are equal i.e. 
7 + =7-, the state will irretrievably lose the coherence, 
but the coherence can be recovered when the magnitudes 




FIG. 2: Time optimal evolutions of a mixed state governed 
by the master equation with the Lindblad operators having 
magnitudes (71,72,73) = (1,0,0) in the descending order of 
the eigenvalues of ((26} . Curves starting from r = (0,0,0.8) 
and approaching r = (0,0, —1) are trajectories of the Bloch 
vector in the x-z plane cross-section of the Bloch sphere. The 
initial direction of the curves is different for each initial angle 
between r and s, i.e. r ■ s = rs cos n = 0, 1, .., 5. 
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As a final remark for this section, we would like to 
point out that, in order to make the time duration of 
the transition physically well-defined, one can introduce 
a small but finite error region around the target state. 
That is, one can be more interested in reaching the tar- 
get state with a fixed fidelity close to one. Then, for 
example, while mathematically it takes an infinite time 
for the system to reach the target state, the system actu- 
ally approaches the target state and enters into its small 
surrounding region in a finite time. We will show in Sec. 
VI that in some cases the time optimal mixed state evo- 
lution can be faster than the time optimal pure state 
evolution (thick solid meridian curve in Fig. 1) discussed 
in [H. 



V. A MODEL FOR MEASUREMENT OR 
DECOHERENCE 

To get a further physical insight into our formalism, 
we study a simple model consisting of two interacting 
spins, one of which is identified with our system (A) and 
the other is an externally controllable ancilla (B). The 
extra ancilla spin can be regarded as representing either 
a measurement apparatus or an environment. The two- 
qubit Hamiltonian can be chosen without loss of gener- 
ality (modulo local rotations of the system and ancilla 
qubits) as 



H AB (t) := ^2 h jk{t)<Jj ® o k , 



(36) 



j,k=x,y,z 



with time-dependent, tunable couplings hjk{t) [301. We 
simulate the optimal Hamiltonian and Lindblad opera- 
tors discussed in the previous section by tuning the cou- 
plings in p6|) and the ancilla state. The state of the 
system is described by the density operator (TT9"]) and the 
state of the ancilla by 



Ps(t) 



1 



[l + b(t)-tr], 



(37) 



where 6 is a tunable Bloch vector. The master equation 
in the Lindblad form |T]) comes from the repetition of 
the unitary evolution by the interaction H AB for a short 
time duration r (i.e., for t much smaller than the typical 
dynamical timescale of p(t)) and the partial trace over 
the B state. Namely, 



P (t + T) = Tr B [e- iHABT p{t) ® p B (t)e iHABT ] 
= p{t)-iTTr B [H AB ,p(t)® PB {t)] 

■< Tr B [H AB (p(t) ® p B {t))H AB ) 

.1 



Tr B [-{H AB ,p{t)® PB {t)}] 



+ 0(t 3 ). 



(38) 



For our model Hamiltonian we can perform the partial 
trace over the B state and get 



p{t + r) = p(t)-ir[H,p(t)} 
+ r 2 

j,k=x,y,z 

+ 0(t 3 ), 



a 3 p(t)o k - ^{(T k (Tj,p(t)} 



(39) 



where the effective single qubit Hamiltonian H is given 
in terms of the couplings hjk and the Bloch vector bj as 



H(t) := ^ hjkbko 

j,k=x,y,z 



J ■ 



(40) 



while the Lindblad matrix is defined by 



a jk{t) ■ = ^ hjihkr, 
l,m—x,y^z 



Sim - 



,—x,y,z 



(41) 

The eigenvectors and the eigenvalues of (|4"Tj) together de- 
fine the components of the Lindblad operators. There- 
fore, since the same components are also the eigenvectors 
in (|26[) . one has to impose the commutativity condition 
[a,K] = 0, where K is given by (|27[) . The real and imagi- 
nary parts of the latter condition give three and nine con- 
straints, respectively, so that we have twelve constraints 
in total for the nine parameters of hjk and the three 
parameters of bj. Then we can adjust the Hamiltonian 
couplings hjk and the components bj of the Bloch vector 
of the ancilla state so that the optimal r(t), s(t), h(t) and 
l a {t) obtained in the previous section are reproduced. 

In particular, if r and s arc parallel (and chosen to be 
along the z-axis as at the end of the previous section), we 
can see that it is sufficient to choose the coupling matrix 
hjk 



(hj k ) 



p 
pOO 
q 



(42) 



with real numbers p, q, and the Bloch vector b as 

b = be 3 . (43) 
The Lindblad matrix (14T1) becomes 



(ajk) 



p 2 ibp 2 
-ibp 2 p 2 



q 2 
Then the K matrix for 6 = 0, 



(44) 



2s 



-i 

1 
r 



(45) 



and the Lindblad matrix a can be diagonalizcd simultane- 
ously. In particular, the eigenvalues of a are p 2 (l T b), q 2 
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and the corresponding eigenvectors are given as in (|32|) 
and (|33p , with the magnitudes of the Lindblad operators 



To 2 



v. 



(46) 
(47) 



The optimal solution for r(t) can be finally cast in the 
form 



r(t) = [-b+{r(0) + b)e-^ 2Tt }e 3 . 



(48) 



An intuitive explanation for the behavior of the Bloch 
vector is the following. The exponential damping in (|48|) 
can be attributed to the standard formula for the tran- 
sition probability calculated from the interaction Hamil- 
tonian 



Hab = p(cr x ®v y J r(Jy® a x ) + qa z ® a z . 



(49) 



Note that q does not appear in (|48|) so that it can be put 
equal to zero. Suppose now that the B spin is up, i.e. 
6 = 1. Since the interaction Hamiltonian (|49|) is propor- 
tional to er+ <g> (7+ — (7_ ® cr_ , both the system and the 
ancilla spins tend to align down. Suppose we start with 
the completely mixed state r(0) = 0. In the statistical 
interpretation of the density operator, this means that in 
half of the systems the spin is up while in the other half 
the spin is down. The up components are steered down 
while the down components remain unchanged, so that 
eventually the state of the system evolves from a com- 
pletely mixed state to a purely down state. This kind 
of 'cooling' has been already discussed before [3l|. On 
the other hand, for b = 0, the system approaches a com- 
pletely mixed state regardless of the initial conditions, as 
we expect because the system is kept into contact with 
the completely random state of the ancilla, ps = 1/2. 
This corresponds to decoherence. We would like to re- 
mark that, although the coupling parameters hjk and 
the ancilla state are constant, in our particular parallel r 
and s case, they can be time dependent in general. This 
particular demonstration of a simple two qubit model suf- 
fices to illustrate the physics behind the master equation 
in the Lindblad form. 



VI. OPTIMAL UNITARY EVOLUTION PLUS 
FINAL MEASUREMENT 

So far we have considered the situation in which the 
evolution of mixed states is the result of a series of short 
time unitary operations each followed by a measurement. 
However, the evolution of mixed states can be also de- 
scribed by a finite time unitary evolution in an enlarged 
Hilbert space for the original system plus an environ- 
ment followed by the trace over the environment, which 
is equivalent to a final measurement only. In this section 
we are going to study in more details this single measure- 
ment strategy and the dependence of the optimal dura- 
tion time of the evolution on the size of the environment. 



To be more explicit, let us consider a set of n-qubits in 
which the first qubit is regarded as the system and the 
rest is the environment. The whole closed system is then 
governed by the Schrodinger equation with a controllable 
Hamiltonian H(t) and the evolution of the first qubit can 
be expressed as 



p(t) = Tr B {U(t)p(p)®pBU*(t)}, 



(50) 



where pB is the initial state of the environment and 
U(t) = Texp(-« / H(t)dt). 

Let us assume that there are no other restrictions for 
the system except the normalization conditions (JTDJ) and 
(fT5|) for the Lindblad and the Hamiltonian operators, re- 
spectively. Now, to justify the dependence of (| L8[) for the 
traceless Hamiltonian on the dimension N (for n qubits, 
N = 2 n ) of the Hilbert space over which it acts, we note 
that 



1 

N 



TrH z 



(51) 



is required by the invariance under the trivial extension 
H^H® I M , i.e. j^- Tr{H ® I M ) 2 = ^ TrH 2 . This is 
because the physical properties of the system described 
by H should not change even if we add an extra ancillary 
environment of arbitrary dimension M which is not inter- 
acting with the system. A simliar argument applies to the 
Lindblad operators. In particular, the constraints should 
be invariant under the trivial extension of the Hilbert 
space. This simple observation is crucial to obtain the 
correct relation between the size of the environment and 
the optimal duration time of the evolution. The same 
point will also become important in addressing the effi- 
ciency of quantum computation and its scaling properties 
when we apply our formalism to explicit quantum algo- 
rithms (32| . 

Therefore, the quantum brachistochrone equation Q16p 
leads to a time independent optimal Hamiltonian for the 
whole closed system which, without any loss of generality, 
can be chosen as 



H = V2^ T cj(|0)(2"- 1 | + ^""^(Ol), 



(52) 



where the states |0) and |2 n_1 ) are the states ||) ® | j) <S> 
... (8> |t) and |J.) ® ||) ® ... ® ||), respectively, in the binary 
representation. The time evolution of the whole system 
is driven by the unitary matrix 



U(t) = cos(V2"- 1 wt)(|0)(0| - 
- ishnV2"-W)(|0)(2 r 



1 >(2- 1 |) 
^XOI). (53) 



Suppose now that we want to steer the first qubit from 
the state p(0) = | T) (T I up to the target state 1 1) (| | ■ Then, 
assuming that p B = (|T)(T|)®""\ from (50]) and |J53J) we 
obtain 



Pit) ^ \ 



1 



cos 



(2v / 2^^)(|T)(T|-|l)ai) 



i sm 



(2V2^i)(|T)UMIXTI) 



,(54) 
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so that the time required to reach the target state is 



T = 



(55) 



which is faster than the pure state case, where the time- 
optimal evolution from |f) to ||) under the same con- 
straint takes the time T = ^- 121. We conclude the 
section by stressing the fact that the optimal duration 
time (|51[) is an exponentially decreasing function of the 
number of the ancilla qubits of the environment. 



VII. SUMMARY AND DISCUSSION 

We have developed a framework based on the vari- 
ational principle for finding the time optimal quantum 
operation to make a transition between given initial and 
final states pi and pf, when the physical system obeys a 
Markovian master equation in the Lindblad form. The 
equations of motion for the Hamiltonian H and the Lind- 
blad operators L a can be written down once the con- 
straints for H and L a are specified according to the 
problem. One then obtains the time optimal operation 
(H(t), L a (t)) and the optimal duration time T by solv- 
ing the quantum brachistochrone and the other equations 
of motion and imposing the initial and final conditions 
p(0) = pi and p(T) = p f . 

Our formalism for the mixed state case has been ex- 
plicitly demonstrated with a one-qubit model. First, 
the optimal Hamiltonian was obtained from the quan- 
tum brachistochrone equation. Then the time evolution 
of the density operator and of the Lindblad operators 
which represent an optimal measurement or decoherence 
was found from a remaining set of ordinary differential 
and algebraic equations. In a particular case an analyti- 
cal solution was given, while the solution for more general 
situations was shown numerically To get a more phys- 
ical intuition, we constructed an interacting two-qubit 
model where an ancilla qubit plays the role of the envi- 
ronment and we demonstrated that repeated short-time 
Markovian transitions can reproduce the optimal time 
evolution of mixed states. Next we considered the time 
optimal evolution of mixed states driven by a single final 
measurement after the unitary evolution in an enlarged 
Hilbert space. In this case the optimal duration time 
is an exponentially decreasing function of the number of 
the qubits of the environment. 

Let us compare the efficiency of the time optimal evo- 
lutions between given initial and final pure states for the 
two models of repeated Markovian measurements and 
that of a final measurement following a unitary evolu- 
tion in an extended Hilbert space. In Fig. 3 we plot 
the fidelity between the target pure state and the time 
optimal evolved mixed state as a function of the dura- 
tion time of the evolution. We can see that, for a given 
fidelity close to one, the evolution via a final measure- 
ment becomes more and more efficient (i.e., it takes a 



' 1 

FIG. 3: Time dependence of the fidelity squared (i|p(i)|J.)- 
The dashed and the solid curves correspond to the repeated 
Markovian measurement case (|48[) and the single, final mea- 
surement case (|54j> in an n-qubit system (n — 1,2,3), re- 
spectively. Until the two curves cross each other, the re- 
peated Markovian measurement model approaches the target 
Pf = 1 1) (I I faster than the single measurement model. 



shorter time) than the evolution via a repetition of mea- 
surements as the number of qubits in the environment 
increases. This is because with an environment more re- 
sources are available for the processing of the information 
required to generate the time optimal evolution. 

Incidentally, we note that both the dynamical evolu- 
tion law (JTJ) (or, equivalcntly, (|39|l ) and (|50|) can be also 
expressed as a completely positive, trace preserving map 

V(t), i.e. p{t) = V(t)p(0) := T,C Wa(t)p(0)WKt), 
where the W a are the Kraus operators, which satisfy 

J2a=o WaW a = 1. In particular, the relation between the 
Lindblad operators in ([T]) and the Kraus operators is ex- 

plicitly given by W = 1 - J J2a=o L l L a and W a = yjrh a 
(see, e.g., (U). 

Our work has not dealt with the more general case 
of different duration times for the contacts between the 
system and the environment and the case of the possi- 
ble memory feed-back from the environment itself. The 
authors of [3l| also considered the problem of control in 
dissipative quantum dynamics in order to achieve optimal 
purification of a quantum state, but they worked within 
the standard framework of a set of constant Lindblad 
operators. Furthermore, although there should be no 
conceptual difficulty in extending our work to the prob- 
lem of optimal quantum control via quantum feedback 
by introducing a stochastic term in the master equation 
[H, HH HI] > we nave not discussed this problem here. 
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